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^ ! Abstract 

By using the fact that Polychronakos-like models can be obtained through the 
'freezing limit' of related spin Calogero models, we calculate the exact spectrum as 

^ • 

well as partition function of SU{m\n) supersymmetric Polychronakos (SP) model. It 
turns out that, similar to the non-supersymmetric case, the spectrum of SU(m\n) SP 
model is also equally spaced. However, the degeneracy factors of corresponding energy 
levels crucially depend on the values of bosonic degrees of freedom (m) and fermionic 
degrees of freedom (n). As a result, the partition functions of SP models are expressed 
through some novel g-polynomials. Finally, by interchanging the bosonic and fermionic 
degrees of freedom, we obtain a duality relation among the partition functions of SP 
models. 
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1 Introduction 



Integrable one- dimensional models with long ranged interactions have recently at- 
tracted a lot of interest due to their close connection with diverse subjects like fractional 
statistics, random matrix theory, level statistics for disordered systems, Yangian alge- 
bra, q-polynomials etc. [1-20]. The SU(M) Polychronakos spin chain [6-8] is a well 
known example of such integrable models with the Hamiltonian given by 

Hp = E J— —^2 > (1-1) 

l<i<j<N [ x i x j) 

where e = 1(— 1) represents the ferromagnetic (anti- ferromagnetic) case and Py is the 
exchange operator interchanging the 'spins' (taking M possible values) of z-th and j- 
th lattice sites. Moreover the positions of corresponding lattice sites (xj), which are 
inhomogeneously distributed on a line, are given by the zero points of the iV-th order 
Hermite polynomial and may also be obtained as a solution of the following set of 
equations: 

* = E 7 _ 2 v» ' (L2) 

k=/=i \ X i X k) 

where i G [1, 2, • • • , N]. Though the Polychronakos spin chain (|1 . 1| ) does not enjoy 



translational invariance, one can find out its exact spectrum as well as partition function 
by considering the 'freezing limit' of the Calogero Hamiltonian which possesses both 
spin and particle degrees of freedom [8]: 

^ i=\ \ ux i I l<i<j<N \ X i X 3) 

u and / being some positive coupling constants. Thus, it is revealed that the Hamilto- 
nian ( |1 . 1| ) generates an equidistant spectrum, where the energy levels are highly degen- 
erate in general. Such high degeneracy of energy levels can be explained through the 
'motif picture originating from the Y(qIm) Yangian symmetry of the Polychronakos 
spin chain ( |1 . 1|) and Calogero model (|1.3|) [13,21]. Moreover, the partition function of 
the Hamiltonian is found to be closely related [13,22] to the SU(M) Rogers-Szego 
(RS) polynomial, which have appeared earlier in different contexts like the theory of 
partitions [23] and the character formula for the Heisenberg XXX spin chain [24]. 



Now, for the purpose of obtaining a supersymmetric extension of the Polychronakos 
spin chain, we consider a set of operators like Cj a (Cj a ), which creates (annihilates) a 
particle of species a on the i-th lattice site. Such creation (annihilation) operators are 
defined to be bosonic if a G [1, 2, • • • , m] and fermionic if a G [m + 1, m + 2, • • • , m + n] 
(where, according to our notation, m + n = M). Next, we focus our attention only to 
a subspace of the related Fock space, for which the total number of particles per site 
is always one: 

M 

E <Afiia = 1, (1-4) 
a=l 

for all i. On the above mentioned subspace, one can define a supersymmetric exchange 
operator as [18] 

M 

f,{m\ri) _ (~,\ n \ n ,-. r\ 



a, /3=1 



and show that this P/ 7 m ' n ' ) yields a realisation of the permutation algebra given by 



V% = 1, VijVji = VuVij = VflVu, [Pij,Vi m ] = 0, (1.6) 

(i, j, I, m being all distinct). This supersymmetric exchange operator ( |1.5| ) was used 
earlier for constructing the SU(m\n) supersymmetric Haldane-Shastry (HS) model [18]. 
So, in analogy with the case of supersymmetric HS model, we may use the exchange 
operator (|1.5|) for constructing a Hamiltonian of the form 

nP n) = £ V tJ .J . (i.7) 

l<i<j<N \ x i x j) 

It should be noted that, at the special case m = M, n = 0, i.e. when all degrees of 



freedom are bosonic, P/ J m ' n ' ) (O)) becomes equivalent to the spin exchange operator 



Pij which appears in the Hamiltonian ( |1.1| ). Therefore, at this pure bosonic case, 
( |1 . 7| ) would reproduce the ferromagnetic Polychronakos spin chain ( |1.1| ) where e = 1. 
Similarly, for the case m = 0, n = M, i.e. , when all degrees of freedom are fermionic, 
( |1 . 71) would reproduce the anti-ferromagnetic Polychronakos spin chain ( |1 . 1|) where 
e = — 1. So, when both bosonic and fermionic degrees of freedom are involved (i.e., 
when both m and n take nonzero values), we may say that 7i^ 1 ' n ' ) (|1.7|) represents the 
Hamiltonian of SU(m\n) supersymmetric Polychronakos (SP) model. 



In this article, our aim is to study the spectrum and partition function of the 
above defined SU(m\n) SP model. To this end, in Sec. 2, we introduce an appropriate 
extension of spin Calogero model ( |1.3|) , which would generate the SP model (|1.7|) at 
the 'freezing limit'. In this section, we also demonstrate that the spectrum of such 
extended spin Calogero model is essentially the same with the spectrum of decoupled 
SU(m\n) supersymmetric harmonic oscillators, where each oscillator has m bosonic as 
well as n fermionic spin degrees of freedom. In this way we are able to show that, the 
equidistant spectrum of SU(m\n) SP model Q1.7Q can be obtained by simply 'modding 
out' the spectrum of SU(m\n) harmonic oscillators through the spectrum of spinless 
bosonic harmonic oscillators. Subsequently, in Sec. 3, we derive the partition function of 
SP model (|1.7| ) and interestingly observe that such partition functions can be expressed 
through some novel g-polynomials. Finally, we obtain a duality relation between the 
partition functions of SU(m\n) and SU(n\m) SP models. Sec. 4 is the concluding 
section. 



2 Spectra of SU{m\n) SP model and related spin 
Calogero model 

For obtaining the spectrum of SP model, we wish to follow here the approach of Ref.8 
and find out at first a suitable extension of spin Calogero model (|1.3j ), which would 
reproduce the SP model (|1.7| ) at the 'freezing limit'. However, one immediate problem 



which arises at present is that a spin Calogero model like ( |1.3|) is a first quantised 
system, while the SP model flL7|) is a second quantised system. So, for applying the 
above mentioned approach, it is convenient to transform the second quantised SP model 



( |1 . 7| ) to a first quantised spin system. 

To this end, we now consider some special cases of 'anyon like' representations of 
permutation algebra ( |1.6| ), which were used earlier for constructing integrable exten- 
sions of SU(M) Calogero-Sutherland model as well as HS spin chain [25-28]. Such 
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special cases of 'anyon like' representations are defined by their action on a spin state 
like \a1a2 ■ ■ ■ un) (with ctj G [1, 2, • • • , M] ) as 

1 ; |aia 2 • • • «i • • • a, • • • on) = e i«»>«»+i> w \a x a2 • • • a,- • • ■ on ■ • • ajv) , 

(2.1) 

w h ere $( m W( a . ai+1) ... )aj ) = Oifai, aj G [1, 2, ■ • • , m], $W n ')(ai,ati +1 ,--',a j ) = 
7r if ati, acj G [m+1, m+2, • • • , m+n], and $ (m|n) (ai, a i+x , ■■■ ,otj) = 7r E™=m+i Ep=*+i 
5 TjQ , p if aj G [1, 2, • • • , m] and <x,- G [m + 1, m + 2, • • • , m + n] or vice versa. It is clear 
that P^ ' reproduces the original spin exchange operator and P^ M ^ reproduces 



-Pij. Next we notice that, due to the constraint ( |1.4| ), the Hilbert space associated 



with SP Hamiltonian ( |1.7| ) can be spanned through the following orthonormal basis 
vectors: C\ ai C2 a2 ■ ■ ■ C* N |0), where |0) is the vacuum state and G [1,2, ■ ■ ■ , M]. 
So, it is possible to define a one-to-one correspondence between the state vectors of the 
above mentioned Hilbert space and the state vectors associated with a spin chain as 

l«i« 2 • • -a N ) <-> C\ ai Cl 2 ■ ■ ■ tf NaN |0) . (2.2) 

However it should be noted that, the matrix elements of P^J 1 ^ ( |2.1|) and P^J 1 (|1.5| ) 
are related as [28] 

(o/^ ■ ■ ■ a' N \P lj \aia 2 ■ ■ ■ a N ) = (0|C 7 v Q / jv C A r_i i ^_ i • ■ • C Xa >P\^ n) C\ ai C\ a2 ■ ■ ■ C j Na jO), 

(2.3) 

where c^s and a[s may be chosen in all possible ways. Thus one finds that, the 'anyon 

like' representation ( |2.1| ) is in fact equivalent to the supersymmetric realisation (|1.5| ). 

Consequently, the first quantised spin Hamiltonian given by 

( 1 _ p( m \ n ) \ 

H? n) = £ 1 - ~- 2 ' , (2-4) 



l<i<j<N 



would also be completely equivalent to the second quantised SP model ( |1.7|) . In par- 
ticular, the Hamiltonians ( |1.7| ) and ( |2.4| ) would share the same spectrum and their 
eigenfunctions can be related through the correspondence ( |2.2|) . 
Next, we define a spin Calogero model as 

1 n / pa \ I ( I - F (m|n) "l 

* 1=1 \ udj t / l<i<j<N \ X i X j) 



where P- ™ is the 'anyon like' representation ( |2 . 1| ) . At the special case m = M, n = 
(m = 0, n — M), the above Hamiltonian reproduces the SU(M) spin Calogero model 
( |1.3|) with e = 1 (e = —1). Notice that the Hamiltonian ( |2.5|) can be rewritten as 

H { ^ n) = Ho + lHi m]n) , (2.6) 

where Hq is the Hamiltonian for spinless Calogero model: 

z i=l \ UJj % J l<i<j<N \ X i X 3) 

and H[ m ^ is obtained from H^ n ^ ( [2^4]) by replacing and with x« and Xj respec- 
tively. Thus, the operator Hi" 1 ^ possesses both spin and particle degrees of freedom. 
However, in analogy with the case of usual spin Calogero model ( |1.3| ) [8] , we may now 
consider the 'freezing limit' of extended spin Calogero model Q2.5p . Such 'freezing limit' 
is obtained by setting I — uj in the Hamiltonian ( |2.5[ ) and finally taking its uj — > oo 
limit. It is evident that, at this 'freezing limit', particle and spin degrees of freedom of 
Hamiltonian (|2.5|) would decouple and the operator H[ m ^ would be transformed to a 
pure spin Hamiltonian where the fixed values of XjS are determined through the minima 
of potential energy associated with spinless Calogero Hamiltonian (|2.7|) . Since the solu- 
tion of eqn.( |1.2|) leads to the minima of such potential energy, the operator H[ m ^ will 
exactly reproduce the SP Hamiltonian ( p.4|) at the 'freezing limit'. Consequently, the 
eigenf unctions of spin Calogero model ( |2.5|) will factorise into eigenfunctions of spin- 
less Calogero model ( |2.7|) containing only particle degrees of freedom and SU(m\n) SP 
model ( |2.4| ) containing only spin degrees of freedom. Moreover, at this 'freezing limit', 
the energy eigenvalues (denoted by E PtS (uj), where the subscripts p and s represent 
the particle and spin degrees of freedom respectively) of the full system ( |2.5[ ) can be 
expressed as 

E p , a (u) = E p (u) + uE 8 , (2.8) 

where E p (uS) and E s denote the energy eigenvalues of spinless Calogero model (p.7| ) 
and SP spin chain ( |2.4|) respectively. 

It is clear from eqn.( |2.8|) that, by 'modding out' the spectrum of spin Calogero 
model ( |2.5|) through the spectrum of spinless Calogero model ( p.7|) , one can construct 
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the spectrum of SP model (|2.4|) or (|1.7|) . So, for the purpose of obtaining the spectrum 



of SP model, it is essential to find out at first the spectrum of spin Calogero model 



Q2.5 ). To this end, we consider the spinless Calogero model of distinguishable particles, 
which is described by a Hamiltonian of the form 

^ i=l \ ^1 / l<i<j<N l x « x jJ 



1 ^ , ' ' ' 1 X 



1 ^Ji 



where is the coordinate exchange operator: Kijip(- 
if)(- ■ ■ , Xj, ■ ■ ■ , Xi, ■ ■ •). By restricting the action of the above Hamiltonian on com- 
pletely symmetric wave functions, we may recover the spinless Calogero model ( |2.7| ) 
of indistinguishable particles. It has been recently found that, by using some simi- 
larity transformations, one can completely decouple all particle degrees of freedom of 
Calogero Hamiltonians ( |2.7| ) as well as ( |2.9| ) [29-33]. Such similarity transformations 
naturally lead to a very efficient method of calculating the eigenf unctions of spinless 
Calogero models. In particular, the spectrum and eigenf unctions of Ho (|2.9| ) can be 
obtained through a similarity transformation which maps this interacting Hamiltonian 
to a system of decoupled harmonic oscillators like 

TV 

H free = u ^ a\a k , ( 2 -10) 

k=l 

where a\ = -7= (u)x k — gf^-J and a k = -7= \0Jx k + The above mentioned simi- 

larity transformation is explicitly given by [32,33] 

S- 1 {Hq -E g ) S = H free , (2.11) 

where E g = \Nu ( Nl + (1 - I) ) and 

S = ^e-^e^ 2 e^ x \ (2.12) 

with 

N N q2 I ^ N 

x2 = Y, x h v2 = E-^2> = II 1%-^l'exp -^E^ 2 

3=1 j=l ux j l<j<k<N \ z j=l 

N d 2 , ^ f 1 ( d d \ K jk - 1 1 



^ L J~[ dx) + 1 ^ { (xj - x k ) \dxj dx k ) + (xj - x k ) 



The operator Ol is called the Lassalle operator. As is well known, the eigenfunc- 
tions of decoupled oscillators (|2.10|) can be written in the form: \n1ri2 ■ ■ ■ TCjv) = 
T\f=i ( a ]) 3 |0), where |0) is the corresponding vacuum state. The coordinate rep- 
resentation of such eigenfunctions, having eigenvalues ujY,jLi n j, is given by 

1 2 N 

ip ni ,n 2 ,-,n N (xi,x 2 , ■ ■ ■ ,x N ) = e~2 wX Y[H nj (Xj), (2.13) 

J'=l 

where H n . (xj) is the Hermite polynomial of order rij. It should be noted that, the above 
eigenfunctions do not generally obey any symmetry property under the exchange of 
coordinates and, therefore, they are free from any statistics. Due to the similarity trans- 
formation ( 2.11 ), the eigenfunctions of spinless Calogero model ( |2.9| ) may be obtained 
through the eigenfunctions ( |2.13| ) of free oscillators as 

Xn U n 2 ,-,n N (x 1 ,x 2 , - ■ ■ ,x N ) = S ^e~^ x2 J| H n .(Xj) \ , (2.14) 

with E{ ni y = E g + uj J2f=i n j representing the corresponding eigenvalues. So, apart 
from a constant energy shift (E g ) for all levels, the spectrum of spinless Calogero model 
( |2.9|) is exactly the same as the spectrum of decoupled Hamiltonian ( 2.1 0|) containing 
distinguishable particles [32,33]. A comment might be in order. The exponentiation of 
the Lassalle operator e~^° L yields essential singularities at = Xj, i,j = l,2,---,N, 
when it operates on general multivariable functions. However, such essential singular- 
ities are not generated when the exponentiation of the Lassalle operator operates on 
multivariable polynomials. Our using of the Lassalle operator in the expression ( |2.14j ) 
is such a safe one. Details on the property of the Lassalle operators are given in refs. 
32, 33. 

In the following, we want to show that the eigenfunctions of spin Calogero model 
(|2.5|) can also be obtained through the eigenfunctions of decoupled harmonic oscillators, 
provided they possess nondynamical spin degrees of freedom and obey some definite 
statistics. To this end, we introduce a projection operator A^ 1 ' r ^ which satisfies the 
relations like 

K^P^A^ = A^KijpW") = A^ |n) , (2.15) 



where i, j G [1, 2, ■ • • , N]. Such a projector can be formally written through 'transpo- 
sition' operator rA m (= P^^Kij) as 

A Hn) = ^ ^ 7 g^... 7 ^ ) (2.16) 

where the series of transposition t^^t^J 1 ^ ■ ■ ■ r^ n ^ represents an element of the per- 
mutation group (-P/v) associated with iV objects and, due to the summations on p 
and {ikijk]i each element of P/v would appear only once in the r.h.s. of the above 
equation. For example, A 2 m ' n ^ and A 3 m ' n ^ are given by: A^""^ = 1 + t[™^ , 

« (mln.) -i . (mln) (mln) . (m\n) (mln) (m\n) (mln) (mln) T , i u u a j 

A3 = 1 + T12 + r 23 + T i3 + r 23 T i2 + r i2 r 23 -It should be noted 
that at the special case m = M, n = (m = 0, n — M), the projector A^ 1 ' n ' 1 ( 2.16Q 
completely symmetrises (antisymmetrises) any wave function under simultaneous in- 
terchange of particle as well as spin degrees of freedom, and thus projects the wave 
function to bosonic (fermionic) subspace. Next, we multiply the eigenfunction ( |2.14|) of 
Calogero Hamiltonian (|2.9j) by an arbitrary spin state \a\a 2 ■ ■ -oln) and subsequently 
apply the projector (|2.16|) for obtaining a wave function like 

K\:Z:"-\Z^) = (7172 •■ -in\ a ( n h U (V^* 2 ft j \^---a N ) 

(2.17) 

where x = x\, x 2 , • • ■ , x/v and 7 = 71, 72, • • • , 7jv- Since the operator Ho ( |2.9| ) commutes 
with Kij as well as A^'™- 1 ( |2.16| ), the expression ( |2.17| ) again gives an eigenfunction 
of the Calogero Hamiltonian ( j2.9|) , where each particle possesses some nondynamical 
spin degrees of freedom. Consequently, with the help of relation Q2.15|) which allows 
one to 'replace' by JF^ 71 , one can show that ^^'nl'.'-'nZi^ 7) ( |2.17|) also gives an 



n 



r 



eigenfunction of spin Calogero model ( |2.5|) with eigenvalue E^ ni y^ ai y = E g +u J2j=i 

Thus we interestingly find that, eqn. (|2.17|) produces all eigenfunctions of spin 
Calogero model ( |2.5| ) through the known eigenfunctions of decoupled harmonic oscil- 
lators ( |2.iq ). However it is important to notice that, due to the existence of projector 
A^'™"* in eqn.( [2.17| ), a definite correlation is now imposed among the eigenfunctions 
of decoupled harmonic oscillators. To demonstrate this point in an explicit way, we 



9 



observe first that the operator S ( 2.12| ) commutes with Kij and the projector A 



(m\n) 
N 



Q2.16|) . By using such commutation relations, eqn.( |2.17|) can be rewritten as 



K\:Z:::.>Z^) = s ni:Z'''Z( x ^)> ( 2 - 18 ) 

where 

KllZPZ&l) = (7i72---7^v|ASv m|n) | (V^ X2 ft Hn^xj)^ \a x a 2 ---a N ) 

(2.19) 

This ^n{]n2l---',riN ( x i l) represents a correlated eigenfunction (with eigenvalue uY^=\ n j) 
for decoupled harmonic oscillators (|2.10|) , where each oscillator carries M number of 
nondynamical spin degrees of freedom. To determine the precise nature of the above 



mentioned correlation, we use the relation ( |2.15| ) and replace A^ n ^ by A ™ KijP^ 



in the r.h.s. of eqn.( [2.19| ). Finally, by acting on \a\0i2 • ■ ■ a/v) and Ky on 



13 x 13 

(m\n) 



e 2 



hujX 2 t-[N 



Ilj=i H Uj .(xj) , we find that the eigenfunction (|2.19| ) of free oscillators must 



satisfy the following symmetry condition under simultaneous interchange of related 
particle and spin quantum numbers: 

^%^:%PZ& 7) = e^ m ^ a ^-^ M/;:;;:;;;; ;: ;: : ;:;;:;;-(.r: 7 ) , (2.20) 

<J>( m l n )(e*i, • • • , a.j) being the same phase factor which appeared in eqn.( [2.1|) . It 
is clear from eqn.( [2.20"D that, for the case ctj, aj G [1, 2, • • • , to], the eigenfunction 
x(x; 7) remains completely unchanged under simultaneous interchange 
of particle and spin quantum numbers. Thus ctj may be treated as a 'bosonic' quantum 
number, if it takes any value ranging from 1 to m. Next, by using eqn.( |2.20] ) for the 
case «j, ay G [m + 1, m + 2, • • • , m + n], it is easy to see that ^nll---lnT,---',n J j',---',nN (x] 7) 
would pick up a minus sign under simultaneous interchange of particle and spin quan- 
tum numbers. Therefore, the eigenfunction ^n^'---'nl,'---,'nj,'---,'n J ^(^; 7) must be trivial if 
we choose = rij and «j = atj G [to + 1, m + 2, • • • , to + n]. Thus «j may be treated 
as a 'fermionic' quantum number, if it takes any value ranging from m + 1 to m + n. 
Notice that if we simultaneously interchange a 'bosonic' quantum number a, with 
a 'fermionic' quantum number atj and rii with rij, then the eigenfunctions satisfying 
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transformation relation fl2.20|) would either remain invariant or pick up a minus sign 
depending on whether even or odd number of fermionic spin quantum numbers are 
present in the configuration: ai + i(Xi + 2 ■ ■ •£*,•_].. It is obvious that, at the special case 
m = M, n = (m = 0, n = M), '"^(x; 7) ( p,19|) represents completely sym- 

metric (antisymmetric) eigenf unctions of SU(M) bosonic (fermionic) oscillators. So, 
for the case m, n 7^ 0, we may say that the correlated state vectors ( |2.19|) would repre- 
sent the eigenf unctions of N number of SU(m\n) supersymmetric harmonic oscillators. 

Due to the existence of symmetry condition (|2.20| ), we find that all independent 
eigenfunctions of SU (m\n) supersymmetric harmonic oscillators can be obtained uniquely 
through the following occupation number representation. Since, at present, spins be- 
have as nondynamical degrees of freedom, all 'single particle states' for this occupation 
number representation may be constructed by taking m + n copies of each energy eigen- 
state associated with a spinless harmonic oscillator — the first m copies being bosonic 
in nature and last n copies being fermionic in nature. As usual, any bosonic single par- 
ticle state can be occupied with arbitrary number of particles and any fermionic single 
particle state can hold at most one particle. By filling up such bosonic and fermionic 
single particle states through N number of particles, we can easily identify all inde- 
pendent eigenfunctions of the form (|2.19|) . So, this occupation number representation 
of states ( p. 19 ) gives us a very convenient way of analysing the spectrum of SU(m\n) 
supersymmetric harmonic oscillators. In particular we may verify that, similar to pure 
bosonic or fermionic case, the spectrum of SU(m\n) supersymmetric harmonic oscil- 
lators is also equally spaced. But, it is interesting to note further that, the ground 
state energy of SU(m\n) supersymmetric harmonic oscillators coincides with that of 
the pure bosonic case, instead of pure fermionic case. Since M number of fermionic sin- 
gle particle states with zero energy can hold at most M number of particles, a nonzero 
ground state energy is obtained for a pure fermionic system when N > M. However, 
if at least one bosonic single particle state with zero energy is available, we can fill up 
that state through all available particles. Consequently we find that, irrespective of 
the values of m and n, the ground state energy of SU(m\n) supersymmetric harmonic 
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oscillators would always be zero. Moreover, by using the above mentioned occupation 
number representation, we obtain the degeneracy (D g ) of ground state for SU(m\n) 
super symmetric harmonic oscillators as 

" (N + m- k-l)\n\ , 
n = i ! i (2 21) 

9 £f (N -k)\(m-l)\k\(n- k)\ ' l ' ' 

It is worth observing that this degeneracy factor crucially depends on the values of 
m and n, and reproduces the degeneracy of the ground state for SU(M) bosonic 
oscillators [8] at the special case m = M, n = 0. Similarly, one can demonstrate 
that the degeneracy of higher energy levels, appearing in the spectrum of SU(m\n) 
supersymmetric harmonic oscillators, would also depend on the values of m and n. 

Due to the relation ( |2.18| ), where S acts as a nonsingular operator, there exists a 
one-to-one correspondence between the independent eigenfunctions of spin Calogero 
model fl2.5|) and SU(m\n) supersymmetric oscillators. Consequently, up to a constant 
energy shift (E g ) for all energy levels, the spectrum of spin Calogero model (2.5) ex- 
actly coincides with the spectrum of SU(m\n) supersymmetric harmonic oscillators. 
However it is well known that, up to the same constant energy shift for all energy 
levels, the spectrum of spinless Calogero model (|2.7|) exactly coincides with the spec- 
trum of iV number of spinless bosonic harmonic oscillators. As the eigenvalues of both 
SU(m\n) supersymmetric oscillators and spinless bosonic oscillators depend linearly 
on the coupling constant uj, it is clear that eqn. fl2.8|) can be used even for any finite 
value of to (though the eigenfunctions of spin Calogero model ( |2.5| ) will factorise only 
at the 'freezing limit'). Moreover, since the spectra of both SU(m\n) supersymmetric 
oscillators and spinless bosonic oscillators are equally spaced, it automatically follows 
from eqn. (|2.8| ) that the spectrum of SU(m\n) SP model ( |1.7|) would also be equally 



spaced for any choice of m and n. But it is natural to expect that, similar to the case 
of SU(m\n) supersymmetric oscillators, the degeneracy of energy levels for SU(m\n) 
SP model (|1.7|) would crucially depend on the values of m and n. So, it should be 
interesting to derive the partition function of SP model (L7|), where all information 
about these degeneracy factors is encoded. 
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3 Partition function of SU(m\n) SP model 

In the previous section we have observed that, similar to the pure bosonic case, the 
ground state energy of SU(m\n) supersymmetric harmonic oscillators would always be 
zero. So, in analogy with this pure bosonic (i.e., ferromagnetic) case [8], we may now 



put u> = 1 in eqn. (|2.8|) and subsequently use this equation to obtain a relation like 

= hJ^-, (3.i) 



where q = e~~^r , Z^ n \q) and Z^ n \q) denote the canonical partition functions of 
SU(m\n) SP model ( |1 . 7|) and SU(m\n) supersymmetric harmonic oscillators (where 
u = 1) respectively. Due to the above mentioned notations, Z^°\q) and Z^\q) 
denote the canonical partition functions of N number of spinless bosonic and fermionic 
harmonic oscillators respectively. It is well known that, the partition functions of such 
spinless bosonic and fermionic oscillators are given by [8] 

(q)n (q)n 

where the standard notation: (q) jy = (1 — q) (1 — q 2 ) ■ • • (1 — q N ) (and (q)o = 1) is used. 
So, for calculating zffl n \q) with the help of eqn.(|3.1|), we have to find out only the 



partition function of iV number of SU(m\n) supersymmetric oscillators. 

To this end, however, we consider at first the grand canonical partition function of 
SU(m\n) supersymmetric oscillators. Such grand canonical partition function may be 
denoted by 2'( m l n ) (q, y), where y = e~ M and \x is the chemical potential of the system. 
So, according to our notation, Z^°\q,y) and Z^ l \q,y) denote the grand canonical 
partition functions of spinless bosonic and fermionic harmonic oscillators respectively. 
As usual, the grand canonical partition function of SU(m\n) supersymmetric oscillators 
can be related to the corresponding partition function through the following power 
series expansion in variable y: 

CO 

Z {mln) (q,y) = E y N ziP ln \q) (3.3) 
JV=0 
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where it is assumed that ZQ m ^ n \q) = 1. In the previous section we have found that, 
all independent eigenf unctions of SU(m\n) super symmetric harmonic oscillators can 
be obtained through an occupation number representation, where the corresponding 
single particle states are constructed by taking m + n copies of each energy eigenstate 
associated with a spinless harmonic oscillator — the first m copies being bosonic in 
nature and last n copies being fermionic in nature. By exploiting this result, it is 
easy to prove that the grand canonical partition function of SU(m\n) supersymmetric 
oscillators can be expressed through those of spinless bosonic and fermionic oscillators 
as 

ZW»\q,y) = [^(Mfp^y)]". (3.4) 

Next, we substitute the power series expansion (|3.3| ) to the place of all grand canonical 
partition functions appearing in the above equation. Subsequently, we compare the 
coefficients of y N from both sides of eqn.( |3.4|) , and readily find that the canonical 
partition function of SU (m\n) supersymmetric harmonic oscillators may also be related 
to those of spinless bosonic and fermionic oscillators: 

m n 

z^ ln \ q ) = e n^%)LK%)' ( 3 - 5 ) 

Y,Zi'«+Y,U b t= N i=1 j=1 

where a^s and bjS are some nonnegative integers. By substituting the known partition 
functions (|3.2| ) of spinless bosonic and fermionic oscillators to the above equation, one 
may now obtain an explicit expression for the partition function of SU(m\n) super- 
symmetric oscillators as 

b o {h j~ 1) 

Z*j ^ ^ (q) = 5 7 — \ — 7 — \ 7 — \ 7 — \ — 7 — \ 7 — \ — ■ (3-6) 

V m a+V" b=N WM(9)«B ,,, (?k , (?)»i(?)i2- ,, (9k 



Finally, by using eqns.( |37LD , ( |3.2| ) and ( p.6| ), we derive the exact canonical partition 
function of SU{m\n) SP model ( |1.7| ) as 

6,(6,-1) 



V m a+V" b=N ' "Worn • Wbi{Q)b 2 ■■■{q)b n 
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Since, there exists an upper bound on the highest energy eigenvalue of SP model (|L 
the partition function Z ™ (q) (T3.7Q evidently yields some new g-polynomials which 
are characterised by the values of m and n. The coefficients of various powers of 
q, appearing in such novel g-polynomial, would represent the degeneracy factors of 



corresponding energy levels associated with SP model ( |1 . 7|) . It is interesting to note 
that, by putting m = M, n = to the expression (|3~7| ), one can exactly reproduce the 
partition function of ferromagnetic Polychronakos spin chain [8]. On the other hand, 
for the limiting case m = 0, n = M, eqn.( [3.7|) would reproduce the partition function 
of anti-ferromagnetic Polychronakos spin chain up to an insignificant multiplicative 
factor. This is because, in pure fermionic system of harmonic oscillators with 

M spin degrees of freedom gives some nonzero ground state energy when N > M. So, 
the r.h.s. of eqn.( p.l| ) must be modified through a multiplicative factor for taking into 
account such nonzero ground state energy [8]. Consequently the expression ( p.7|) for 
partition function, which we have derived for the supersymmetric case (i.e., when both 
m and n are nonzero), should also be modified by the same factor at the pure fermionic 
limit. 

With the help of a relation: (q~ 1 ) l = (— 1)^ g (q)i, we find that the partition 
function (|3.7|) of SP model satisfies a remarkable duality condition given by 



ZiP ln \q) = q^ ZP m \q^) , (3.8) 

where m and n may be chosen as any nonzero integer. Due to this duality condition, 
one can write down a relation of the form 

D< m|n) (£) = DP m) ( ^^f X) ~ , (3.9) 



where DpJ 1 (E) denotes the degeneracy factor associated with energy eigenvalue E of 
SU(m\n) SP model. Moreover, by using the duality condition (|3.8|) , along with the 
fact that the ground state energy of Hamiltonian ( |L7|) is zero, it is easy to show that 



Bm „ = EKpH ( 3,o) 
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represents the highest energy eigenvalue of the SU(m\n) SP model. In this context one 
may notice that, the highest energy eigenvalue of ferromagnetic or anti-ferromagnetic 
SU(M) Polychronakos spin chain is given by [8] 

M - 1 „ t(M -t) . 
E max = N - -i '- , 3.11 

where t = N mod M. Thus we curiously find that, in contrast to the case of SU(M) 
Polychronakos spin chain (|1 . 1|) , the highest energy eigenvalue of SU(m\n) SP model 
( |1.7| ) does not depend at all on the values of m or n. 

For obtaining some insight about the above mentioned difference between the high- 
est energy eigenvalues ( |3.10| ) and (|3.11 ), we finally consider the 'motif picture, which 
was used to analyse the degeneracy of eigenfunctions for SU(M) HS as well as Poly- 
chronakos spin chain through the corresponding symmetry algebra [11,13]. As is well 
known, motifs of these spin chains are made of binary digits like '0' and l V. Moreover, 
for a spin chain with A" number of lattice sites, these binary digits would form motifs 
of length N — 1. So we can write a motif as (a±a2 ■ ■ ■ fljv-i), where Oj G [0, 1] and each 
motif represents a class of degenerate eigenfunctions which yield an irreducible repre- 
sentation of Y(glM) Yangian algebra [11,18]. For the case of SU(M) Polychronakos 
spin chain, the energy eigenvalue corresponding to (aia 2 ■ ■ ■ o-n-i) motif is given by [13] 

N-l 

^Kaa-a^.!) = T a r ' ( 3 ' 12 ) 

r=l 

However, for the case of SU(M) Polychronakos and HS spin chain, there exists a 
'selection rule' which forbids the occurrence of M number of consecutive 'l's in any 
motif. By combining this selection rule along with eqn. fl3.12|) , one can derive the 
highest energy eigenvalue ( |3.11[ ) of SU(M) Polychronakos model and also understand 
why this eigenvalue crucially depends on the value of M. It is worth noting that, the 
above mentioned motif picture can also be used to analyse the spectrum of SU(m\n) 
super symmetric HS spin chain [18,28]. But, for this supersymmetric case, there exists 
no 'selection rule' and binary digits like '0' and '1' can be chosen freely for constructing a 
motif of length N—l. Being motivated by the case of supersymmetric HS spin chain, we 
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may now conjecture that the motifs corresponding to SP model (|1.7j) are also free from 
any 'selection rule' and eqn. (|3.12|) again gives the corresponding energy eigenvalues. 
So, the spectrum of SU(m\n) SP model contains 2 N ~ 1 number of motifs, which can 
be constructed by filling up iV — 1 positions with '0' or '1' in all possible ways. By 
using eqn.( p.l2|) , it is easy to check that the motif (11 • ■ • 1) yields the highest energy 
eigenvalue (|3.10| ). Thus the absence of any 'selection rule', for the motifs corresponding 



to SP model, turns out to be the main reason behind the remarkably simple expression 

(ELD- 



4 Concluding Remarks 

In this paper we have investigated the spectrum as well as partition function of SU (m\n) 
supersymmetric Polychronakos (SP) model ( |1.7| ) and the related spin Calogero model 
( |2.5| ). The similarity transformation ( |2.11| ), which maps the spinless Calogero model 



of distinguishable particles to decoupled oscillators, and the projection operator ( |2.16| ) 
have played a key role in our derivation for the spectrum of spin Calogero model fl2.5p. 



Thus we have found that, up to a constant energy shift for all states, the spectrum of 
this spin Calogero model is exactly the same with the spectrum of decoupled SU(m\n) 
supersymmetric harmonic oscillators. Furthermore, by using the occupation number 
representation associated with SU(m\n) supersymmetric harmonic oscillators, we have 
obtained the exact partition function for this spin Calogero model. 

It turned out that, the above mentioned spin Calogero model reproduces the SP 



model (L7) at the 'freezing limit'. Consequently, by factoring out the contributions 
due to dynamical degrees of freedom from the spectrum as well as partition function 
of spin Calogero model Q2.5|), one can compute the spectrum and partition function for 
SU(m\n) SP model. By following this procedure we have found that, similar to the 
non-supersymmetric case, the spectrum of SU(m\n) SP model is also equally spaced. 
However, the degeneracy factors of corresponding energy levels crucially depend on the 
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values of m and n. As a result, we get some novel (/-polynomials which represent the 
partition functions of SP models. Moreover, by interchanging the bosonic and fermionic 
degrees of freedom, we obtain a duality relation among the partition functions of SP 
models. 

As a future study, it should be interesting to find out the Lax operators and con- 
served quantities for the SU(m\n) SP model (|1.7| ). Moreover, in parallel to the case 
of SU(m\n) Haldane-Shastry model [18], one might be able to show that the SU(m\n) 
SP model also exhibits the Y^gl^n)) super- Yangian symmetry. Such super- Yangian 
symmetry of SP model may turn out to be very helpful in analysing the degeneracy 
patterns for its spectrum. However, information about these degeneracy patterns is 
also encoded in our expression of partition function (|3.7|). So, there should exist an 
intriguing connection between the partition function ( ^.7[ ) and the motif representa- 
tions for super- Yangian algebra. In particular, the partition function ( |3.71) might lead 
to a supersymmetric generalisation of well-known Rogers-Szego (RS) polynomial. The 
recursion relation among these supersymmetric RS polynomials may then be used to 
find out the motifs representations and degenerate multiplets associated with super- 
Yangian algebra. We hope to report about such supersymmetric RS polynomials and 
related motifs in a forthcoming publication [34]. 
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